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[57] ABSTRACT 

The present invention provides a turbulence model that has 
been developed for blade-element helicopter simulation. 
This model uses an innovative temporal and geometrical 
distribution algorithm that preserves the statistical charac- 
teristics of the turbulence spectra over the rotor disc, while 
providing velocity components in real time to each of five 
blade-element stations along each of four blades, for a total 
of twenty blade-element stations. The simulator system 
includes a software implementation of flight dynamics that 
adheres to the guidelines for turbulence set forth in military 
specifications. One of the features of the present simulator 
system is that it applies simulated turbulence to the rotor 
blades of the helicopter, rather than to its center of gravity. 
The simulator system accurately models the rotor penetra- 
tion into a gust field. It includes time correlation between the 
front and rear of the main rotor, as well as between the side 
forces felt at the center of gravity and at the tail rotor. It also 
includes features for added realism, such as patchy turbu- 
lence and vertical gusts in to which the rotor disc penetrates. 
These features are realized by a unique real time implemen- 
tation of the turbulence filters. The new simulator system 
uses two arrays one on either side of the main rotor to record 
the turbulence field and to produce time -correlation from the 
front to the rear of the rotor disc. The use of Gaussian 
Interpolation between the two arrays maintains the statistical 
properties of the turbulence across the rotor disc. 

The present simulator system and method may be used in 
future and existing real-time helicopter simulations with 
minimal increase in computational workload. 

20 Claims, 8 Drawing Sheets 
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FIGURE 3B 
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SYSTEM AND METHOD FOR FINITE 
ELEMENT SIMULATION OF HELICOPTER 
TURBULENCE 

ORIGIN OF THE INVENTION 

The invention described herein was made in the perfor- 
mance of work under a NASA contract and is subject to the 
provision of Section 305 of the National Aeronautics and 
Space Act of 1958, Public Law 85568 (72 Stat. 435; 42 
U.S.C. 2457). 

BACKGROUND OF THE INVENTION 

1. Technical Field of the Invention 

The present invention is generally directed at aircraft 
simulation systems. More specifically, this invention relates 
to a new system and method for creating turbulence, includ- 
ing progressive gust penetration by the rotor, for use in 
real-time piloted helicopter simulations. 

2. Description of the Prior Art 

The modern helicopter is no longer a vehicle used for 
simple missions where only its hovering capability is 
required. The helicopter, and rotorcraft in general, are 
increasingly faced with complex missions which push the 
aircraft to its design limits. One emerging requirement 
which challenges helicopters is stabilized flight through 
moderate or severe atmospheric turbulence to accomplish 
high workload mission tasks. 

Rotorcraft, both civilian and military, now compete with 
many other forms of transportation and have often shown 
great reliability and productivity in mission performance. 
These rotorcraft may be required to perform missions such 
as hovering over, and landing on moving platforms in high 
turbulence conditions. In light of the dynamic flight condi- 
tions which are encountered routinely, it is essential to have 
the capability to accurately predict rotorcraft performance in 
a turbulent atmosphere. 

It is therefore clear that there is a great and still unsatisfied 
need for real-time rotorcraft atmospheric turbulence models 
in flight simulators. One of the most common grievances 
that many pilots have is that while they do feel turbulence 
acting on the helicopter in the flight simulator, the motion is 
not similar to real atmospheric turbulence. 

One of the most severe limitations of the current rotorcraft 
atmospheric turbulence models is the lack of understanding 
of the fundamental physics underlying the rotorcraft 
response to atmospheric turbulence, which is different from 
that of the fixed wing aircraft. There is a fundamental 
difference between the ways in which rotary wing aircraft, 
such as helicopters, and fixed wing aircraft, such as 
airplanes, experience atmospheric turbulence. The differ- 
ence is primarily due to blade rotational velocity. Whereas a 
wing element of an airplane has pure translational motion as 
the airplane cuts through a turbulence field, a blade element 
of a helicopter has translational as well as rotational motion. 
In an airplane both the turbulence and the response to 
turbulence are stationary random processes which require 
conventional body-fixed sampling, By comparison, the tur- 
bulence as well as the response to turbulence in a helicopter 
are cyclostationary and require blade-fixed sampling. 

The rotational velocity introduces appreciable spatial dis- 
tribution of turbulence velocities over the rotor disc. For 
conventional helicopters operating at attitudes of 1000 feet 
or below, the rotational velocity effects are not negligible. 
Research efforts have been made to provide methods with an 
analytical basis to predict and simulate turbulence and 
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response statistics. One such simulation method is based on 
representing turbulence sample functions in terms of multi- 
variable sinusoids with random phases and provides second- 
order statistics of covariances and spectral densities. 

5 Another turbulence simulation method is described in 
Costello, M. F., “A Theory for the Analysis of Rotorcraft 
Operating in Atmospheric Turbulence,” Proceedings of the 
46th Annual National Forum of the American Helicopter 
Society, Washington, D. C., May 1990, pp. 1003-1015. This 
10 article describes turbulence over a rotor blade that is 
approximated by a series of radial shape functions. With this 
expansion, a stochastic state space model is formed where 
the system dynamics matrix and the control matrix are 
periodic with a period equal to the rotor rotational speed and 
15 input to the plant is generated by independent white noise 
sources. 

The following articles reflect some of attempted methods 
for improving turbulence simulation models, all of which 
are incorporated herein by reference: 

Howlett, J. J., “UH-60A Black Hawk Engineering Simu- 
lation Program: Vol. I -Mathematical Model,” NASA CR-1 
66309, Dec. 1981. 

Prasad, J. V. R., Riaz, J., Gaonkar, G. H., Yingyi, D., “Real 
25 Time Implementation Aspects of a Rotorcraft Turbulence 
Simulation Method,” 49th Annual Forum of the American 
Helicopter Society, St. Louis, Mo., May 1993, p. 459. 

Dahl, H. J. and Faulkner, A. J., “Helicopter Simulation in 
Atmospheric Turbulence,” Vertica, Vol.3, 1979, pp. 65-78. 
30 Judd, M. and Newman, S. J., “An Analysis of Helicopter 
Rotor Response due to Gust and Turbulence,” Vertica, Vol. 
1, 1977, pp. 179188. 

Dang, Y., G. Gaonkar, Prasad, J., Zhang, H., “Parallel 
Computing of Helicopter Response to Turbulence Toward 
35 Real-Time Implementation,” 50th Annual Forum of the 
American Helicopter Society, Washington, D.C., Vol. 2, 
May 1994, pp. 869-882. 

However, the foregoing conventional methods have not 
proven to be totally satisfactory in real-time simulations 
since they either take excessive execution time or they do 
not yield realistic turbulence according to pilot opinion. 

SUMMARY OF THE INVENTION 

45 The present invention provides a unique simulator system 
and method for use in future and existing real-time helicop- 
ter simulations with minimal increase in computational 
workload. 

The present invention further provides a turbulence model 
50 developed for blade-element helicopter simulation. This 
model uses an innovative temporal and geometrical distri- 
bution algorithm that preserves the statistical characteristics 
of the turbulence spectra over the rotor disc, while providing 
velocity components in real time to each of five blade- 
55 element stations along each of four blades, for a total of 
twenty blade-element stations. 

An initial investigation of the simulator system was 
performed using a piloted, motion-based simulation of the 
Sikorsky UH60A Black Hawk helicopter. Although only the 
60 vertical component of stochastic turbulence was used in this 
initial investigation, vertical turbulence components induced 
vehicle responses in all translational and rotational degrees 
of freedom of the helicopter. 

The single-degree-of-freedom configuration of the 
65 present simulator system was compared to a conventional 
full 6-degrees-of-freedom baseline configuration, where 
translational velocity inputs are superimposed at the vehicle 
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center of gravity, and rotational velocity inputs are created 
from filters that approximate the immersion rate into the 
turbulent field. 

The simulator system includes a software implementation 
of flight dynamics that adheres to the guidelines for turbu- 
lence set forth in military specifications. One of the features 
of the present simulator system is that it applies simulated 
turbulence to the rotor blades of the helicopter, rather than 
to its center of gravity. The simulator system accurately 
models the rotor penetration into a gust field. It includes time 
correlation between the front and rear of the main rotor, as 
well as between the side forces felt at the center of gravity 
and at the tail rotor. It also includes features for added 
realism, such as patchy turbulence and vertical gusts into 
which the rotor disc penetrates. 

Briefly, the foregoing and other previously unattainable 
features and advantages of the present invention are realized 
by a unique real time implementation of the turbulence 
filters. The new simulator system uses two arrays, one on 
either side of the main rotor, to record the turbulence field 
and to produce time -correlation from the front to the rear of 
the rotor disc. The use of “Gaussian Interpolation” between 
the two arrays maintains the statistical properties of the 
turbulence across the rotor disc. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other features of the present invention and 
the manner of attaining them will become apparent, and the 
invention itself will be best understood, by reference to the 
following description and the accompanying drawings, 
wherein: 

FIG. 1 is an isometric view of a helicopter cockpit module 
for use in conjunction with the simulator system and method 
according to the present invention; 

FIG. 2 is a top view of a simulated rotor disc forming part 
of the simulator system and method of the present invention, 
showing the rotor velocities and the added turbulence 
velocities; 

FIGS. 3A, 3B, 3C and 3D represent a functional block 
diagram of the present simulator system titled “TURBR”; 

FIG. 4 represents a block diagram of the Patches of 
Turbulence circuit used in the present simulator system; and 

FIG. 5 represents a block diagram of the Vertical Gusts of 
Turbulence circuit used in the present simulator system. 

Similar numerals refer to similar elements in the drawing. 
It should be understood that the sizes of the different 
components in the Figures are not in exact proportion, and 
are shown for visual clarity and for the purpose of expla- 
nation. 

DETAILED DESCRIPTION OF THE 
INVENTION 

As used herein, turbulent fluid motion is an irregular 
condition of flow in which the various quantities show a 
random variation with time and space coordinates, so that 
statistical distinct average values can be discerned. Accord- 
ing to this definition, the flow has to satisfy the condition of 
irregularity which highlights an important feature of turbu- 
lence. Because of the irregular nature of any turbulent flow, 
it is very difficult to describe it by deterministic functional 
relationships. On the other hand, turbulence velocities are 
generally random in the sense that it is possible to describe 
it by the laws of probability, hence a statistical approach is 
often adopted to describe atmospheric turbulence velocities. 

The flight simulator used for performing the piloted 
simulation according to the present invention is NASA 
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Ames Research Center’s Vertical Motion Simulator (VMS). 
This simulator is a general purpose research tool that retains 
a high degree of flexibility to enable tailoring for a wide 
range of fixed and rotary wing applications. The simulator 
5 can be configured to meet the needs of a particular task by 
selecting hardware and software options. In a specific 
experiment related to the present invention, the simulator 
was configured for the purpose of evaluating the influence of 
stochastic turbulence created at the blade-element stations of 
10 a helicopter rotor simulator system. 

Some key facilities of the simulator system are comprised 
of a motion simulator, a visual simulator, pilot controls, 
primary flight instrumentation, and a software or processor 
implementation of flight dynamics. These components are 
15 briefly described below. 

The motion simulator system generally provides a cockpit 
module 10 with motion in all six axes (roll, pitch, yaw, 
heave, lateral and longitudinal), as shown in FIG. 1. The 
motion simulator system is capable of high accelerations, 
20 velocities and displacements, and notably, the performance 
described herein can be achieved simultaneously in all six 
degrees of freedom. 

The visual simulator system generally provides a visual 
representation of the flight environment through monitors 
25 mounted in the cockpit module 10. The cockpit module 10, 
pilot controls and primary flight instrumentation are similar 
to those of a GENHEL UH60A Black Hawk helicopter. The 
software implementation of flight dynamics including all the 
requisite motion and visual communication software, has 
been executed on a single processor, with a cycle time of 
about 12 milliseconds. 

The experimental simulations conducted at NASA Ames 
Research Center met with very favorable pilot opinion. 
35 Some of these simulations were implemented into the GEN- 
HEL UH60A mathematical model described in Howlett, J. 
J., “UH-60A Black Hawk Engineering Simulation Program: 
Vol. 1- Mathematical Model,” NASA CR-166309, Dec. 
1981. An heuristic model of turbulence was developed that 
40 preserves the Gaussian statistical characteristics of turbu- 
lence filters over the rotor disc. This model avoids the large 
computational and storage requirements of conventional 
investigations into rotorcraft turbulence models, such as 
those described in Prasad, J. V. R., Riaz, J., Gaonkar, G. H., 
45 Yingyi, D., “Real Time Implementation Aspects of a Rotor- 
craft Turbulence Simulation Method,” 49th Annual Forum 
of the American Helicopter Society, St. Louis, Mo., May 
1993, p. 459. 

The present simulator system is designed to create and 
50 propagate all three translational components of turbulence in 
real time. The turbulence velocity components are superim- 
posed at each blade element location as a function of the 
helicopter velocity, with distribution according to the blade 
element instantaneous geometry. 

55 The entire stochastic rotor element turbulence model will 
be described herein. Examples from this model are presented 
for the helicopter in low speed flight (i.e., 10 knots). In this 
study the Dryden spectral form is selected due to its com- 
putational efficiency, although a discrete model of a curve-fit 
60 approximation to the von Karman form may alternatively be 
used. Spectral analysis of the simulation shows that the 
rotor-to-body attenuation is so severe that the original gust 
spectra appear to be of secondary consideration. 

Conventional military specifications (“MIL SPEC”) 
65 angular turbulence filters are a substitute for a finite-element 
distribution of translational turbulence in the air mass over 
the airframe. In the present simulator system only the rotor 
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disc may be considered as an “airframe”, and the transla- 
tional turbulence in the air mass is distributed over the rotor 
disc to each blade element. This distribution produces rotor 
moments and forces. The MIL SPEC translational turbu- 
lence filters may be used in the present simulator model. 5 

Referring to FIG. 2, it is a diagrammatic representation of 
the rotor geometry, illustrating a rotor disc 12 in top plan 
view. The helicopter horizontal velocity component in the 
forward direction is represented by V H , and the added 
turbulence velocity components are represented by Au L , 
Av l , Aw l , Aur, Av^, Aw^. Since the rotor disc 12 encounters 
turbulence at its leading edge 14, the line of turbulence 
which translates along the rotor disc 12 based on the velocity 
of the helicopter, is represented by the onset line 16. The 
onset line 16 is perpendicular to the vehicle forward flight 1 
path (i.e., the direction of V^,), and tangential to the leading 
edge of the rotor disc 12. The vertical added turbulence 
velocity components Aw z and Aw^ and the horizontal added 
turbulence velocity components Au z , Av i? AUr, and Avj, 
which will be described later in greater detail, are rep re- 20 
sented by corresponding arrows at the left and right ends of 
the onset line 16. 

From Taylor’s hypothesis, the turbulent flow is stationary 
and homogeneous in the longitudinal direction. Turbulence 25 
filter outputs are statistically valid at either fixed spatial 
locations or at points translating with the helicopter, as is the 
case in the conventional body-fixed formulation. The out- 
puts of the turbulence filters are applicable at points on the 
onset line 16 at the discrete time points. Interpolation 30 
perpendicular to the onset line 16 requires a transport delay 
that is proportional to the distance along the flight path to the 
blade element location, and inversely proportional to the 
helicopter aerodynamic velocity. The turbulence is assumed 
uncorrelated at each side of the rotor disc 12, and homoge- 35 
neous in the lateral dimension. The statistical properties of 
the turbulence are then invariant over the entire rotor disc 12 
by using “Gaussian Interpolation” between values appli- 
cable on each side of the rotor disc 12. 

40 

In a particular simulation using the Black Hawk helicop- 
ter model, the rotor system has four blades (N=4). For each 
of the blades the number of segments used in the simulation 
was five (M=5). The blade and segment indices are defined, 
as follows: 


l<K<iV (blades) 

1 <m<M (segments) 

The finite element analysis is performed using in part the 
equal-annuli algorithm described in Howlett, J. J., “UH-60A 
Black Hawk Engineering Simulation Program: Vol. I 
-Mathematical Model,” NASA CR-1 66309, Dec. 1981, to 55 
generate the radii (r m ) to the blade stations. r m may be 
computed as follows: 



where Q=27 rad/sec is the nominal rotor RPM. 

In the present model, rather than defining the translational 
disturbance velocities at the helicopter center of gravity, they 
are defined at two onset points 20, 21 located along the onset 
line 16. Using Gaussian inputs, three independent turbulence 
velocity histories are computed at each of these two points 
20, 21 in real time. These histories are updated each cycle 
time, and the pertinent values from these histories are 
distributed to the individual segment locations by temporal 
and geometrical algorithms. 

The helicopter horizontal velocity component in the for- 
ward direction (v^) may be defined approximately by the 
following equation, where u b is the helicopter longitudinal 
aerodynamic velocity and v b is the lateral aerodynamic 
velocity: 


V H = N Ub 2 + Vb 2 

Using the velocity v H , the sideslip angle “p” and the 
perpendicular onset line 16 shown in FIG. 2, the aerody- 
namic azimuth angle (wind axis) for the n th blade would 
then be defined by the following equation: 

w n =w n =$ 

where x P n is the n th blade geometrical azimuth angle from 
aft. 

This geometry is sufficient to establish the longitudinal 
distribution. The left and right Dry den filters receive Gaus- 
sian inputs, so the outputs of these linear filters are also 
Gaussian. Furthermore, a transport delay in these outputs 
does not alter the statistical characteristics. The outputs of 
the filters are stored in tables, which are updated each cycle 
time. The value applicable for a given element is determined 
by selecting a value from the tables. Given that the table size 
is established (i.e., K^=500) for storing time histories of the 
outputs of the filters, the minimum aerodynamic velocity 
V min in that may be accommodated is determined by dis- 
tributing the table over the rotor diameter: 

2 R 

V/nin ~ , . . 

K M &t 

If a smaller minimum velocity is required, a larger table 
size must be selected. For the simulation cycle time 
At=0.012 sec, the minimum velocity for K^=500 is 8.944 
ft/sec (5.3 knots). The rotor-plane aerodynamic velocity 
used in the turbulence filters is restricted to this minimum. 

J Vmin Vff < v min 

Vuv= < 

I V H V H = V min 


( m -|) 

r m =\ (e + e') 2 + — [R 2 ~(e + e’) 2 ] - e 

where the rotor radius is R=268 ft, the hinge offset is e=1.2 
ft, and the spar length is e'=2.2 ft. 

The azimuth angle of the n th blade with respect to the 
helicopter aft centerline is given by the following equation: 


For this simulation, independent of vehicle velocity, 
60 exactly K JV/ At=6 seconds of data is stored in the tables. Time 
histories of the outputs of six turbulence filters are contained 
in six separate tables, each of length K M . At the minimum 
velocity the individual cells correspond to a length of 
0.10732 ft (such that 500 cells of the tables span the rotor 
65 diameter), whereas at a velocity of 100 ft/sec (59.2 knots), 
for example, the cells are each 1.2 ft long such that only 45 
cells are needed to span the rotor diameter. 



5,860,807 


8 


7 


In order to select the velocities to be used from the tables, 
for any blade element, an integer index in the table for the 
blade element turbulence (k m? n ) is computed to determine 
the blade element perpendicular distance from the onset line. 
Defining, 5 


A ps = 



7/6 

N 0.8v w 



) 


r m ,„=R+r m cos W n 

the delay index (k m „) is defined by the following equation: 



where the upper brackets denote the integer ceiling opera- 
tion. The tables extend beyond the rotor disc for higher 
vehicle velocities than the minimum, such that stationary 
outputs are available for table size increases caused by 
decreases in velocity. 

The Dry den model will now be described in more detail. 
It is known that the low-altitude vertical turbulence scale 
length (L w ) is given as a piecewise continuous function of 
altitude (h), as follows: 


L w = 


< 


10 h < 10 ft 
h 10 ^ h ^ 1000 ft 
1000 h > 1000 ft 


and the in-plane scale lengths L u and L v are given by the 
following functions: 


L U =L V = 


75.64 h < 10 ft 

hfj, 2 10 ^ h ^ 1000 ft 

1000 h > 1000 ft 


where f uv is defined by: 
/„ v =0. 177+0. 000823/* 


A?0) = 


4b 


■ sAw(s) 


10 


3tV MV 

4b 


Ar(s) = 


3b 


■ sAv(s) 


%v uv 

3b 


15 


These angular disturbances are not used in the illustrated 
simulation model because the geometry is considered in the 
points of application of the translational turbulence excita- 
tions. The asymmetric turbulence velocities produce 
20 moments, which, in turn, produce angular activity. The 
angular disturbence equations are meant to approximate the 
body dynamics. The new turbulence implementation uses 
the body’s own interactions with the rotor to produce 
angular disturbances. 

25 Continuous transfer functions involving random inputs 
are typically converted to discrete form using the zero-order 
hold formulation, where the input is assumed constant over 
each computer cycle At. The z-transforms of Laplace func- 
tions f-(s) then become, 

30 

Fi(z)=z { (^— )/;c, } 


which may then be converted to difference equations for 
35 discrete implementation. 

Using this technique on the Dryden transfer functions 
produces the following coefficients: 


40 


Y u = v„ v A t/L u 
fi = e_Yu 


The horizontal turbulence RMS intensities o u ,a v are a 
function of the selected vertical turbulence intensity a w . 


/2 = a„(l-/i)\T2/r7 




-0.4 


45 


Yv = v„ v Af/L v 


The results shown herein are for o w =5 ft/sec, h=200 ft, 
and v^l.6.9 ft/sec (10 knots). 

Using Taylor’s frozen field hypothesis, the Dryden form 
of the linear turbulence transfer functions are given by the 
following equations: 


Am ( 5 ) = ■ 


oN 


2v w /nL u 


- V U y/L U 


gj = 2e-^ 

= -e- 2Yv 

g3 = [<J V / M"yv][ 1 - e-Yv + (NT7 - l)e-YvY v ] 

g4 = - [a v e-Yv/\[Y^][l - e-Yv + (NTT - 1)y v ] 


Yw = v uv AtjL w 


Av(^) = 


O'y^l Vyy/j XiLy (M7 5 + Vuv/Ly) 

(s + VjLy ) 2 


hj = 

h 2 = - e - 2Yw 


Aw(s) : 


C^wNl V u y/jlLt w " (NTT 5 + V u y/L w ) 

(S + Vuy/Lw) 2 


The angular transfer functions are developed from partial 
differential equations relating these translational distur- 
bances to the immersion rate of the aircraft into the turbulent 
field. These equations use a parameter “b”, which is the wing 
span length of the helicopter. In an exemplary simulation “b” 
is approximately 26 ft. 


60 


h$ = [<J W / M Yw ][1 - e~Yw + (NTT - l)e _ YwY w ] 
/*4 = — [a^e-Y^/Nl Yw ][1 - e'YH- + (NTT - l)Yw] 


These coefficients are used in the following six difference 
65 equations, each driven by an independent Gaussian noise 
source %, with a zero mean value and a unity standard 
deviation. 
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Au l {K)=f t Au L (k-l^n^k) 

Au R (k)=f 1 Au R (k- l)+f 2 n 2 (k) 

5 

Av L (k)=g 1 Av L (k- l)+g 2 Av L (k-2) 

+g 3 n 3 (k)+g 4 n 3 (k-l) 

10 

Av R (k)= gl Av R (k-l)+g 2 Av R (k-2) 

+g 3 n 4 (k)+g 4 n 4 (k-l) 

Aw l (k)=h 1 Aw L (k- l)+h 2 Aw L (k-2) 

+h 3 n 5 (k)+h 4 n 5 (k- 1) 

AW^k^h^W^k-^+h^W^k-l) 20 

+h 3 n 6 (k)+h 4 n 6 (k-l) 

These computationally efficient equations produce stable 25 
outputs for all stable inputs. For the body-fixed formulation 
these velocities are superimposed at the center of gravity. 
However, in the present simulator system, these velocities 
are created at the two onset points 20, 21 and then distributed 
and superimposed at the various blade element locations. 30 
They produce forces and moments because they modify the 
angle of attack of each blade element. 

The foregoing analytical study will now be correlated 
more specifically to the block diagrams in FIGS. 3 A, 3B, 3C, 35 
with reference to the following two tables of variables. 


TABLE I 


Name 

TURBR VARIABLES 
Description 

BETA 

Sideslip Angle 

DISP 

Dispersion - standard deviation of vertical 
turbulence 

DT2 

Cycle Time 

FIRSTELM 

Value of single input for all segments 

FTR 

Distance from CG to TR in feet 

GAINU 

Scale factor for U-body turbulence 

GAINV 

Scale factor for V-body turbulence 

HCG 

Height 

IFIRSTELM 

Use single input for all segments 

INDX 

Index for segment arrays 

ITRBGUST 

Use vertical gusts 

ITRBPAT 

Use Patchy Turbulence 

ITURBRW 

Allow segment inputs 

ITWBD 

Create vertical body turbulence 

IVCG 

Index for CG for V-array 

IVTR 

Index for Tail Rotor for V-array 

Nx 

Number of cells in the main rotor arrays 

PSIR(M) 

Angle to blade #1 in radians 

R, RMR 

Main Rotor Radius 

SEGMID(N) 

Midpoints of the segments (in fractions or 
RMR) 

UB 

X-velocity (ft/s) 

UPGMR(M,N) 

Segment turbulence components 

UTURB 

U - Body Turbulence 

VB 

Y-velocity (ft/s) 

VTURB 

V - Body Turbulence 

VYGTR 

Tail rotor Y-gust component 

WTURB 

W - Body Turbulence 


TABLE II 


EXTRA SECTIONS 

('VERTICAL GUSTS, PATCHES OF TURBULENCE) 

Gust-Name 

Patch-Name 

Description 

TIME, A(303) 

TIME, A(303) 

Total time since beginning operate 

DT2, A(168) 

DT2, A(168) 

Cycle time 

GTLGUS 

PTLSIG 

Time since last change of output 

GTRNDW 

PTRNDW 

Time to wait until next change of output 

GLMBD 

PLMBD 

Average time to wait between changes 

GLMBDMX 

PLMBDMX 

Max Average time 

GLMBDMN 

PLMBDMN 

Min Average time 

GRMPT 

PRMPT 

Time to ramp in the change 

GRMPTMX 

PRMPTMX 

Max Ramp time 

GRMPTMN 

PRMPTMN 

Min Ramp time 

GVMIN 

PVMIN 

Minimum velocity 

GVMAX 

PVMAX 

Maximum velocity 

GUSTNEW 

PSIGNEW 

New output value to ramp toward 

GOUT 

POUT 

Current output value 


Given the height above ground (h) and vertical standard 
deviation (a w ), the remaining turbulence specifications can 
be defined as follows: The function £ h appendix is given by 
the following equation: 

/ a =0.177+0.000823/j 

TABLE III 


TURBULENCE SPECIFICATIONS 




[ 37.82 

o 

VII 



I hf h ~^/2 

10 < h ^ 1000 

L U ,L V 

s 

h/2 

1000 < h ^ 175 



875 

h ^ 1750 



* 5 

h ^ 10 

Lw 


h/2 

10 < h ^ 175 



. 875 

h > 1750 


1 

f (0.18523)a w 

o 

VII 

a u , a v 

i 

aJV 0 - 4 

10 < h < 1000 


1 

l a w 

h ^ 1000 


FIG. 3 A shows 11 parameters that are input to the 
simulator system, and generally grouped into four different 
groups, as illustrated in the following table IV: 

TABLE IV 


PARAMETERS GROUPING 


Group Parameter 

Rotor Geometry SEGMID(N), PSIR(M), R, FTR, BETA 

Helicopter Dynamics UB, VB, HCG 

Simulator System Variable DT2 

Turbulence Variable Determined NX, DISP 

by User 


FIGS. 3 A, 3B, 3C and 3D combined represent a func- 
tional block diagram of the present simulator system herein 
titled “TURBR”. For simplicity of description, the circuits 
shown in FIGS. 3 A and 3B will be collectively referred to 
as the input circuit 100, the circuit shown in FIG. 3C will be 
referred to as the blade-element turbulence circuit 102, and 
the circuit shown in FIG. 3D will be referred to as the 
bodyturbulence circuit 104. 

Referring to FIGS. 3A and 3B, the input circuit 100 
generates the following output signals that are fed to the 
blade-element turbulence circuit 102 in FIG. 3C: 

AA: corresponds to DT2. 

AB: corresponds to V w , which is calculated as follows: 

v hv =(UB 2 +VB 2 ) 112 
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AC: corresponds to SEGMID(N). 

AD: corresponds to PSIR(M). 

AE: is calculated as follows: 

tan -\VBIUB) 

AF: corresponds to V DTI , which is calculated as follows: 
Vv DT , where v dt—(PT2)-(Y hv) 

AG: corresponds to R. 

The following output signals in FIG. 3B are fed to the 
body -turbulence circuit 104 in FIG. 3D: 

AH: corresponds to and is fed to the Patches of 
Turbulence circuit 150 (FIG. 4). 

AJ: corresponds to Lw in the above Table IV - Turbulence 
Specifications. 

AK: corresponds to Lu in the above Table IV - Turbulence 
Specifications. 

AL: corresponds either to a w at the output of the Patches 
of Turbulence circuit 150 in FIG. 4, or to DISP. 
ITRBPAT is a flag which chooses the “Patches of 25 
Turbulence” a w or a user specified field a w in the 
variable display. 

AM: corresponds to a u in the above Table IV - Turbulence 
Specifications. 

30 

AN: corresponds to a v in the above Table IV - Turbulence 
Specifications. 

AO: corresponds to Lv in the above Table IV - Turbulence 
Specifications. 

AP: corresponds to IVCG, which is calculated as follows: 35 
IVCG=(Nv-l)x0.5, 
where Nv=[(NxlR)(FTR )]+ 2 

AQ: corresponds to IVTR, which is calculated as follows: 40 
[FTRxv on xcos(BETA)]+IVCG. 

A loop 109 performs the following routine: If V^^ is 
greater than, or equal to 2R, no recompute is needed. If, on 45 
the other hand, V DTN is less than 2R, the recompute routine 
110 is executed. This corresponds to the velocity being 
below the minimum velocity allowed by the table size, and 
performs the corresponding equation above relating to Vuv. 

With reference to FIG. 3C the input signals AA, AB, AC, 50 
AD, AE, AF, AG to the blade-element turbulence circuit 102 
are processed as follows: 

AA and AB are input to a vertical gust circuit 200 (shown 
in detail in FIG. 5) whose output may used when a flag 
(ITRBGUST=1), or is not used when the flag 55 
(ITRBGUST=0), and which corresponds to the newest 
of TRBGUST array. 

AD and AE are added to provide the angle of each blade 
(one blade at a time) relative to the direction of translation 
W H . The cosine of the blade angle is calculated and multi- 60 
plied by r m (where r m =ACxAG, as described above) and also 
by AF, to generate INDX. INDX is then simultaneously fed 
to three arrays W A , W R , W z , where the array W A holds the 
outputs of the Vertical Gusts circuit 200, and the arrays W R 
and W z hold the turbulence values found at the onset point. 65 
The sine of the blade angle is also calculated and is pro- 
cessed as follows to generate P m n : 
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P m ,n=[ {[(Sine of blade angle)*(r ra )FJAG]}x0.5]+0.5. 

Input signals CB and CC are respectively fed from the 
body-turbulence circuit 104 (FIG. 3D) to the arrays and 
5 p and the outputs of arrays and W L are fed to a 
Gaussian Interpolation circuit W B whose output OW 5 is 
defined by the following equation: 

Pm,nWR + (1 ~Pm,n)WL 

10 OW B ■■ r ■■ 

N Pmv + (1 -Pmji) 2 

The outputs of array W A and circuit W B are added to 
15 provide the total turbulence input UPGMR(M,N) to the 
blade segment. It should be noted that the circuit inside the 
rectangle in dashed lines constitutes a loop which is repeated 
for each blade-segment combination (i.e., twenty times in 
the present example). 

With reference to FIG. 3D, the input signals AH, AJ, AK, 
AL, AM, AM, AO, AP, AQ to the body- turbulence circuit 
104 are processed as follows: 

AH, AJ, and a Gaussian random number Z(a w ) are fed 
into a Dryden turbulence filter W c for generating signal CC 
according to the following equation: 

N v uv /2jiL w (NT s + v uv /2L w ) 

(s + Vuv /2L w y 

AH, AJ, and the Gaussian random number z(a w ) are fed 
into a Dryden turbulence filter for generating signal CB 
according to the following equation: 

\| vJ2%L w (NJT s + v uv /2L w ) 

(5 + VuvITLw) 2 

Signals AL, AM and AN are fed into a Gaussian Random 
Numbers calculation filter W H , one of whose outputs z(a M ) 
is fed into a Dryden turbulence filter W £ along with signals 
AH and AK, for generating a signal OW £ according to the 
following equation: 

\ 2v wv /jtL w 

( s + Vuv/Lu) 

Signal OW £ is then multiplied by GAINU to generate output 
signal UTURB. 

Another output z (a v ) of the Gaussian Random Numbers 
calculation filter W^ is fed simultaneously into two arrays 
W^ and W M . Signal AP is the index into array W^, and 
signal AQ is the index into array W^. The output of array 
W^ is fed into a Dryden turbulence filter W^ along with 
signals AH and AO to generate a signal OW £ according to 
the following equation: 

NTTTtZT (NJT S + V U y/2Ly) 

(s + VuvIlLy ) 2 

Similarly, the output of array W M is fed into a Dryden 
turbulence filter W G along with signals AH and AO to 
generate a signal OW G according to the following equation: 
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Nl v uv l'2jtL v (VT S + V U y/2Ly) 

(S + VjTLyf 

5 

Each of the signals OW £ and OW G is then multiplied by 
GAINV, to generate output signals VTURB and VTGTR, 
respectively. 

The two arrays W K and W M are useful for correlating the 
side gusts entering the center of gravity and the side gusts 10 
entering the tail rotor. For instance, if the simulation heli- 
copter were flying forward through a sudden side gust, the 
center of gravity of the helicopter would respond before the 
tail rotor enters the gust, thus providing a realistic response. 

In another situation where the helicopter were flying back- 15 
ward through the sudden side gust, the tail rotor will respond 
to the gust before the center of gravity does. Therefore, the 
two arrays and may delay the response of either the 
center of gravity or the tail rotor relative to each other, based 
on the velocity of the helicopter to provide a realistic 2 o 
response to gust conditions. 

In addition, the present model renders a realistic turbu- 
lence simulation and provides a correlation to the terrain and 
other considerations from the physical environment sur- 
rounding the helicopter. To this end, the present model is 2 s 
provided with multiple spatial inputs, such as at the two 
onset points 20, 21, and at two side input locations: the 
center of gravity and the tail rotor. Turbulence is initiated for 
all the segments at the two onset points 20, 21. VTURB 
signal represents the turbulence velocity that is applied to 30 
the center of gravity, and VYGTR represents the lateral 
velocity in the Y direction that is applied to the tail rotor. 

Referring now to FIG. 4, it represents a block diagram of 
the Patches of Turbulence circuit 150. The waiting time Tw 
is calculated according to the following equation: 35 

Tw=Time-PTLSIG 

A software or firmware routine embedded in the circuit 
150 inquires whether Tw is greater than the random wait 40 
time PTRNDW. If Tw is greater than the random wait time 
PTRNDW, the software runs a routine 152 followed by a 
routine 154. If, on the other hand, Tw does not exceed the 
random wait time PTRNDW, then the software runs routine 
154. 45 

The routine 152 will now be described. Vuv is fed into a 
look up table for generating the corresponding PLMBD, 
which is then fed into another look-up table for generating 
output signal PRMPT. Output signals PTRNDW and 
PSIGNEW are also generated according to the following 50 
equations: 

PTRNDW=PLMBDxIn[0.85(U(0,l))+0.1] 

PSIGNEW=abs(z(DISP* 1.2533)) 55 

where U is a uniform random number, and z is a Gaussian 
random number. 

The routine 154 will now be described. The output signal 
a w is determined according to the following equation: 6Q 

o w =[(POUT-PSIGNEW) (DT2)/(PRMPT)]+PO UT 

Referring now to FIG. 5, it represents a block diagram of 
the Vertical Gusts of Turbulence circuit 200. The circuit 200 
is similar to the Patches of Turbulence circuit 150 (FIG. 4), 65 
with the exception that it includes two routines 202 and 204. 
The routine 202 is similar to the routine 152 (FIG. 4) except 
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that GSIG is used according to the following equation, to 
generate an output signal GUSTNEW: 

GUSTNEW=z(GSIG), 

where z is a Gaussian random number. 

While specific embodiments of the system and method for 
finite element simulation of helicopter turbulence have been 
illustrated and described in accordance with the present 
invention, modifications and changes of the apparatus, 
parameters, use and operation will become apparent to those 
skilled in the art, without departing from the scope of the 
invention. 

What is claimed is: 

1. A simulator system for use in real-time piloted heli- 
copter simulation, the helicopter having a center of gravity, 
rotor disc, a tail rotor, a plurality of blades, each of which 
including a plurality of blade-element stations, the system 
comprising in combination: 

a processor for applying finite element turbulence simu- 
lation to the rotor blades of the helicopter, for creating 
and propagating three translational components of tur- 
bulence in real time, and for superposing turbulence 
velocity components at each blade element station as a 
function of helicopter velocity, with distribution 
according to each blade element instantaneous geom- 
etry; 

an onset line representing a line of turbulence along the 
rotor disc, said onset line being tangential to the rotor 
disc, with said turbulence velocity components being 
applied at two points along said onset line; 

said processor periodically updating said turbulence 
velocity components and distributing the updated tur- 
bulence velocity components to the individual blade - 
segment stations; and 

at least one array W K , W M for correlating gusts entering 
the center of gravity and gusts entering the tail rotor, 
and for delaying the response of either the center of 
gravity or the tail rotor relative to each other, based on 
the velocity of the helicopter to provide a realistic 
response to gust conditions. 

2. The simulator system according to claim 1, wherein 
said processor includes an input circuit, a blade-element 
turbulence circuit and a body-turbulence circuit. 

3. The simulator system according to claim 2, wherein 
said input circuit generates a plurality of output signals (AA, 
AB, AC, AD, AE, AF, AG) that are fed to a body-turbulence 
circuit, and that generally correspond to, and are derived 
from the helicopter rotor geometry and dynamics 
parameters, simulator system variables, and turbulence vari- 
ables. 

4. The simulator system according to claim 3, wherein 
said input circuit further generates a plurality of output 
signals (AH, AJ, AK, AL, AM, AN, AO, AP, AQ) that are fed 
to said blade-element turbulence circuit, and that generally 
correspond to, and are derived from the helicopter rotor 
geometry and dynamics parameters, simulator system 
variables, and turbulence variables. 

5. The simulator system according to claim 4, wherein 
said signals AA and AB are input to a vertical gust circuit 
whose output may used when a flag ITRBGUST=1, or is not 
used when said flag ITRBGUST =0. 

6. The simulator system according to claim 5, wherein 
said signals AD and AE are added to provide a blade angle 
of each blade relative to a direction of translation of the 
helicopter; 

wherein the cosine of said blade angle is calculated and 
multiplied by r m (where r m =ACxAG) and also by AF, 
to generate array index INDX; and 
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wherein INDX is simultaneously fed to an array (W A ) that 
holds the outputs of said Vertical Gusts circuit, and to 
two arrays (W^, W z ,) that hold the turbulence values at 
said onset line. 

7. The simulator system according to claim 6, wherein the 5 
sine of said blade angle is calculated and processed as 
follows to generate P : 

P m /J =[{[(Sine of blade angle) *(r m )]+[AG]}x0.5]+0.5. 

8. The simulator system according to claim 7, wherein 10 
signals CB and CC are respectively fed from said body- 
turbulence circuit to said arrays and W L ; and 

wherein P m n and the outputs of said arrays and W L 
are fed to a Gaussian Interpolation circuit whose output 15 
OW 5 is defined by the following equation: 

Pm,nWR + (1 ~p m ,n)W L 

ow B = 

’ Pm,n + (1 -pmji) 2 20 

9. The simulator system according to claim 8, wherein the 
outputs of said array W A and said circuit W B are added to 
provide a total turbulence input to said blade -segment. 

10. The simulator system according to claim 9, wherein 25 
said total turbulence input is determined for each one of said 
plurality of blade -segments. 

11. The simulator system according to claim 10, wherein 
said signal AH corresponds to a velocity (V^) and is fed to 

a Patches of Turbulence circuit; and 30 

wherein signal AL may correspond either to the output of 
said a Patches of Turbulence circuit or to an input signal 
to said Patches of Turbulence circuit. 

12. The simulator system according to claim 11, wherein 
said Patches of Turbulence circuit calculates a waiting time 35 
Tw according to the following equation: 

Tw=Time-PTLSIG, 

where PTLSIC represents the last time; 40 

wherein if Tw is greater than a random wait time 
PTRNDW, said Patches of Turbulence circuit runs a 
first routine followed by a second routine, but if Tw 
does not exceed said random wait time PTRNDW, then 
said Patches of Turbulence circuit runs said second 45 
routine. 

13. The simulator system according to claim 12, wherein 
a rotor-plane aerodynamic velocity Vuv is fed into a look up 
table for generating a corresponding time delay PLMBD, 
which is then fed into another look-up table for generating 50 
an output signal ramp up time PRMPT. 

14. The simulator system according to claim 13, wherein 

said first routine further generates two output signals random 
time to wait (PTRNDW) and turbulence magnitude 
(PSIGNEW) according to the following equations: 55 

PTRND W=PLMBDxln [0.85 (£7(0, l))+0. 1] 

PSIGNEW=abs (z(DISP * 1.2533)) 

where U is a uniform random number, and z is a Gaussian 60 
random number. 

15. The simulator system according to claim 14, wherein 
said second routine generates an output signal a w that 
represents a turbulence RMS intensity. 

16. The simulator system according to claim 15, wherein 65 
said Vertical Gusts of Turbulence circuit calculates said 
ating time Tw, and if Tw is found to be greater than said 
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random wait time PTRNDW, said Vertical Gusts of Turbu- 
lence circuit runs a first Vertical Gusts of Turbulence routine 
followed by a second Vertical Gusts of Turbulence routine, 
but if Tw does not exceed said random wait time PTRNDW, 
then said Vertical Gusts of Turbulence circuit runs said 
second Vertical Gusts of Turbulence routine; 

wherein a rotor-plane aerodynamic velocity Vuv is fed 
into a look up table for generating a corresponding 
delay time PLMBD, which is then fed into another 
look-up table for generating an output signal ramp-up 
time PRMPT; 

wherein said first Vertical Gusts of Turbulence routine 
further generates two output signals time to wait 
(GTRNDW) and gust magnitude (GUSTNEW) accord- 
ing to the following equations: 

Gri?A®W=PLMB£>xIn[0.85(£/(0,l))+0.1] 

GUSTNEW=z(GSIG), 

where z is a Gaussian random number, and GSIG is standard 
derivation of turbulence; and 

wherein said second Vertical Gusts of Turbulence routine 
generates an output signal o w that represents a turbu- 
lence RMS intensity. 

17. A method for use in real-time piloted helicopter 
simulation, the helicopter having a center of gravity, a rotor 
disc, a rotor tail, and a plurality of blades, each of which 
including a plurality of blade -element stations, the method 
comprising the steps of: 

applying finite element turbulence simulation to the rotor 
blades of the helicopter creating and propagating three 
translational components of turbulence in real time, and 
superposing turbulence velocity components at each 
blade element station as a function of helicopter 
velocity, with distribution according to each blade 
element instantaneous geometry; 

defining an onset line for representing a line of turbulence 
along the rotor disc, said onset line being tangential to 
the rotor disc; 

applying said turbulence velocity components at two 
points along said onset line; 

periodically updating said turbulence velocity compo- 
nents and distributing the updated turbulence velocity 
components to the individual blade-segment stations; 
and 

using at least one array W K , W M to correlate gusts 
entering the center of gravity and gusts entering the tail 
rotor, and for delaying the response of either the center 
of gravity or the tail rotor relative to each other, based 
on the velocity of the helicopter to provide a realistic 
response to gust conditions. 

18. The method according to claim 17, wherein an input 
circuit generates a plurality of output signals (AA, AB, AC, 
AD, AE, AF, AG) that are fed to a body-turbulence circuit, 
and that generally correspond to, and are derived from the 
helicopter rotor geometry and dynamics parameters, simu- 
lator system variables, and turbulence variables. 

19. The method according to claim 18, wherein said input 
circuit further generates a plurality of output signals (AH, 
AJ, AK, AL, AM, AN, AO, AP, AQ) that are fed to said 
blade -element turbulence circuit, and that generally corre- 
spond to, and are derived from the helicopter rotor geometry 
and dynamics parameters, simulator system variables, and 
turbulence variables. 

20. A software implemented real-time piloted helicopter 
simulation, wherein the helicopter has a center of gravity, a 
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rotor disc, a rotor tail, and a plurality of blades, each of 
which including a plurality of blade-element stations, the 
simulation comprising in combination: 

a routine for applying finite element turbulence simulation 
to the rotor blades of the helicopter, for creating and 5 
propagating three translational components of turbu- 
lence in real time, and for superposing turbulence 
velocity components at each blade element station as a 
function of helicopter velocity, with distribution 
according to each blade element instantaneous geom- 10 
etry; 

a simulated onset line representing a line of turbulence 
along the rotor disc, said onset line being tangential to 
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the rotor disc, with said turbulence velocity compo- 
nents being applied at two points along said onset line; 
said routine periodically updating said turbulence velocity 
components and distributing the updated turbulence 
velocity components to the individual blade-segment 
stations; and 

a routine for using at least one array W^, to correlate 
gusts entering the center of gravity and gusts entering 
the tail rotor, and for delaying the response of either the 
center of gravity or the tail rotor relative to each other, 
based on the velocity of the helicopter to provide a 
realistic response to gust conditions. 
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